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As light localization becomes increasingly pronounced in photonic systems with 
less order, we investigate optically induced two-dimensional Fibonacci structures 
which are supposed to be amongst the most ordered realizations of deterministic 
aperiodic patterns. For the generation of corresponding refractive index structures, 
we implement a recently developed incremental induction method using nondiffract¬ 
ing Bessel beams as waveguide formation entities. Even though Fibonacci struc¬ 
tures present slightly reduced order, we show that transverse light transport is 
significantly hampered here in comparison with periodic lattices that account for 
discrete diffraction. Our experimental findings are supported by numerical simula¬ 
tions that additionally illustrate a development of transverse light localization for 
increasing propagation distance. 


1 Introduction 

Order is one of the central properties to characterize complexly structured systems of any 
dimensionality. There is certainly no doubt that periodic systems hold highest order and thus 
periodicity - meaning the invariance in translation for integer multiples of a set of lattice 
vectors - are examined thoroughly for centuries providing nowadays a deep theoretical under¬ 
standing of waves in periodic media. Consequently, universal theories such as the Floquet-Bloch 
theorem were introduced where the model of band structures, for instance, is one prominent 
achievement recovered in numerous disciplines. One of the best established model transfers 
from electron-wave propagation in solids is covered by the Floquet-Bloch theorem to the field 
of photonic lattices with light (or electro-magnetic waves in general) as the propagating wave 
in a periodically modulated refractive index material |1,2[. The existence of band gaps where 
propagation is forbidden as, per definition, no eigenstates exist within these gaps is closely 
connected to a band-structure system. 

Quasi-periodic structures obviously are of less order than periodic structures as they lack in 
short-range order |3,|4]. However, the spectra of both, periodic and quasi-periodic structures 
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only present discrete contributions and are distinguishable by their rotation symmetry at the 
utmost which is limited to 2-, 3-, 4-, or 6-fold for periodic lattices. For quasi-periodic structures 
even higher rotation symmetry can emerge and, more interestingly, band-structure properties 
can be assigned to these systems as well, de facto offering complete band gaps [|J6|. 

Numerous experimental techniques to achieve refractive index modulations such as direct 
laser writing [7| and photo lithography techniques [8] were suggested during the last decades [2]. 
In photorefractive media, for instance, the refractive index is usually modulated by illumina¬ 
tion with structured light |9[. In particular, optical induction of elongated two-dimensional 
(2d) photonic structures can be achieved using so-called nondiffracting beams 1101 where the 
intensity is modulated transversely while being constant in the direction of propagation 1111. 
This technique is highly dynamic as the refractive index modulation achieved with low to 
moderate intensities is reversible and additionally introduces a nonlinear response allowing for 
the realization of discrete soliton formations in periodic photonic structures 112,131. Corre¬ 
sponding to the induction of periodic photonic lattices, the use of quasi-periodic nondiffracting 
writing beams such as a 5-fold Penrose intensity-configuration allows for the optical generation 
of according photonic quasi-crystals 114]. 

In general, nondiffracting beams cover an enormous variety of intensity modulations, ranging 
with decreasing order from periodic 15 to quasi-periodic |11 to random structures 116,171. 


Yet, discrete structures without rotation symmetry are barely feasible with a single-beam in¬ 
duction configuration. We thus presented recently that optical induction techniques can be 
extended to resemble any aperiodic structure by sets of zero-order Bessel beams as waveguide 
formation entities |18j. Applying this incremental induction technique |19,20|, we realized full 
aperiodic so-called Vogel lattices that are prominent examples for structures inspired by nature 
as similar arrangements can be found by plant growth following golden-angle phyllotaxis |2Tj. 
Vogel spirals do not show any periodicity in most cases, though they are created deterministi¬ 
cally as a construction rule assembles the structure iteratively up to an arbitrary quantity of 
elements. 

Of course, there are uncountable examples of deterministic aperiodic structures 1221 holding 
different degrees of order [23] 


Yet, in terms of order the Lebesgue’s decomposition theorem 

Spectra of 


was suggested to categorize structures along their spectral appearance |24,25 


aperiodic structures with highest order including Fibonacci patterns |26,27 show only discrete 
contributions while low-order structures such as Rudin-Shapiro patterns feature continuous 
quasi-white-noise spectra 1221. Accordingly, structures with intermediate order have both, 
singular and continuous spectral properties. The most famous structure among them certainly 
arises from the Thue-Morse sequence |23|. 

The question of how an altering degree of order influences wave (or light) propagation remains 
untouched within this categorization scheme. As localization in the sense of reduced light 
transport becomes more significant in disordered photonic systems |28-|30|, the occurrence of 
localized modes due to a deterministic aperiodicity is ascribed to flat bands along with band 
gaps that unclose with decreased order [5]. In consequence, finding signatures of localization 
already in singular-spectral structures would directly refer to an expected reduced degree of 
order in comparison to periodic lattices with likewise singular spectra. 

In this contribution, we examine waveguide patterns of pure-point spectra by exploring light 
localization effects |3T|32 in optically induced refractive index Fibonacci structures. We benefit 
from a highly dynamic and reversible induction scheme at realizing largely elongated photonic 
structures that offer excellent interaction distances. In Fibonacci arrangements, diverse local 
configurations of lattice sites exist |27|. Hence, we consider an output average for waveguide 
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excitation at different input positions. By comparing our results with a periodic photonic lattice 
configuration where discrete diffraction with a high rate of transport is expected to be found [7j, 
we underline that order is significantly diminished already in aperiodic structures with pure- 
point spectra, disabling light transport and causing power-law weighted localization |33|. 

The manuscript is structured as follows. First, we present our approach of applying one¬ 
dimensional Fibonacci words to create 2d lattices with Fibonacci tiling. After we have described 
the applied experimental techniques and specific parameters that we used for the optical in¬ 
duction, we present our results of single-waveguide excitation and compare these experimental 
data with numerical simulations. Finally, we cross-check our results against light propagation 
events in a regular lattice before a conclusion is drawn in the last section. 

2 Design and induction of 2d Fibonacci lattices 

Our approach to design an aperiodic Fibonacci pattern is to vary the distances of adjacent 
waveguides encoded as Fibonacci words 1271. These sequences are binary, and according to 
the Fibonacci series the n th word is generated by combining the (n — l) st and the (n — 2) nd 
word, such as S n = {S n -iS n - 2 }- Giving the first two words determines the complete set of all 
words. We define two different distances A and B — A/(p where ip = (1 + V^)/2 is the golden 
ratio, and set So = A and Si = AB such that the first five Fibonacci words read as 

<So — A, Si = AB , S 2 — ABA , 

<S 3 = ABAAB , <S 4 = ABAABABA. 

Picking two sub-words of length N from a very long Fibonacci word S n with n^> N starting 
at arbitrary but different elements for two transverse directions, we receive a deterministic 
aperiodic, non-symmetric structure with N x N sites as depicted in Fig. [fja). Here, the typical 
character of a Fibonacci word is present in both orthogonal directions. That is, distance A 
occurs with a probability of 0.62, thus more frequently than B which yields typical structure 
groups: quad, double and single waveguide elements. 

The Fourier transform of the Fibonacci structure as shown in Fig. [ija) gives the spatial 
spectrum whose absolute value distribution S(k x ,k y ) is presented in Fig. [ijb). The singular 
character becomes apparent here but especially in Fig. [jje) where the orange plot represents 
the cross section distribution S(k t ) along the arrow drawn in in Fig. |T](b) . For comparison, the 
corresponding plot for a largely extended Fibonacci grating and thus with extensive resolution 
is given by the black curve. Notice that particular frequency peaks have mutual distance 
relations equal to the golden ratio ip. Frequencies with distinct amplitude peaks can be found 
at k — ip 7X1 k where k is the frequency of the prominent peak settled between tv/ A and tv /B and 
m is a natural number. Thus, in addition to the nomenclature aperiodic, a classification of the 
presented spectrum as singular or pure-point is appropriate, as well 1 25] . 

After we have determined the waveguide positioning scheme resembling a 2d Fibonacci 
arrangement, we apply the incoherent-Bessel beam induction method to preserve the non¬ 
diffracting character of the effective intensity during the writing process. This method implies 
to write every single waveguide with an appropriate nondiffracting Bessel beam of zeroth order 
and structural size g. Determining the structural size g of the Bessel beam intensity settles 
the waveguide diameter which scales proportionally with g. In our considerations, we fix g 
to 13 pm and the effective waveguide distance to = 32 pm. This yields A — 37.5 pm and 
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Figure 1 . (a) Fibonacci lattice with Gaussian beam sites (underlying Fibonacci words indicated along¬ 
side). (b) Spatial spectrum S(k^^k y ) according to lattice shown in (a), (c) Measured effective intensity 
with Bessel beam lattice sites taken by multiple-shot illumination at back face of crystal. Bottom right 
quadrant shows according numerical simulation, (d) Experimental output for plane-wave probing (con¬ 
tours indicate waveguide positions), (e) Plot of (orange) spectrum S(k t ) along direction denoted in (b) 
and of (black) ideal spectrum for extended aperiodic Fibonacci lattice with ^-function lattice sites. 


B = 23.2 pm. A simulation of the resulting transverse effective intensity distribution is given 
in the lower right quadrant of Fig. [ljc) . 

In general, our experimental setup incorporating a set of spatial light modulators (SLM) is 
appropriate to generate any kind of nondiffracting beam 111]. The respective setup scheme is 
presented in Fig. [2] and the induction process corresponds to descriptions given in Ref. |18j. 
We use SLMs to experimentally realize numerically calculated light fields in a particular image 
plane defined by an optical imaging system with a demagnification factor of roughly 1/6. 
Particularly, the PSLM (Holoeye Pluto ) is positioned in real space (related to the image plane) 
in order to modulate incoming plane waves. The entire held information of the desired beam 
is encoded in elaborate diffraction gratings displayed by the PSLM. An ASLM (Holoeye LC-R 
2500) is placed in Fourier space for spectral low-pass filtering reasons. For ah experiments, we 
use a frequency doubled Nd:YAG continuous-wave laser source at 532 nm wavelength. 

To resemble the desired 9x9 Fibonacci structure as presented in Fig. |l](a) , we basically 
change the position of each Bessel writing beam and put its central maximum to the defined 
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Figure 2. Experimental setup for induction of Fibonacci lattices with incoherent Bessel beams. 
(A/P)SLM: (amplitude/phase) spatial light modulator, Cam: CCD camera, A/2: half-wave plate, 
LED: background illumination, (R)P: (rotatable) polarizer, SBN: strontium barium niobate crystal; 
lenses and pin hole not labeled. 


position resulting from the given structure design. A full set of diffraction gratings containing 
the necessary writing light fields is generated by convolution of a sparse matrix (with the entire 
position information) and a basic Bessel beam field distribution. 

For an actual experimental induction it is sufficient to send every Bessel beam diffraction 
grating sequentially to the SLM each generating a corresponding light intensity in the volume 
of interest where a phot or efr active strontium barium niobate (SBN) crystal is placed. During 
the illumination, an external field of 2kVcm -1 is applied to the crystal |9|. This sequential 
induction scheme implies an effectively incoherent superposition of all writing light fields of 
one set. The obvious motivation is that a coherent overlap of contributing Bessel beams would 
cause undesired intensity modulations since off-diagonal terms of the resulting field would be 
nonzero. However, by introducing an effective intensity / e ff, we aim to implement the sum of 
all intensities / e ff = l^k| 2 rather than the absolute square of all fields Ek^k| 2 - 

Simulated and experimental effective intensities are given in Fig. [jjc). The intensity distri¬ 
butions of writing and probing light fields are recorded by a CCD camera imaging system. To 
receive the experimental picture, the output intensities of each writing beam are taken indi¬ 
vidually and added up afterward. This overall intensity pattern displays the effective intensity 
that optically induces the 2d photonic structure. A LED is placed above the crystal to actively 
erase the inscribed structure for further light potential inductions. 

Figure [ijd) shows the output intensity when probing the induced structure with a plane 
wave. Notice, that waveguide positions of quad and double elements cannot be resolved prop¬ 
erly here as one intensity envelope covers clusters of adjacent waveguides for predominantly 
perpendicular probe beam incidence. However, we will see later on that single waveguide ex¬ 
citation indicates accurate induction of waveguide clusters according to the effective intensity 
distribution. 

3 Light propagation characteristics in Fibonacci lattice 

After the induction, we experimentally investigate the influence of the Fibonacci lattice with 
its diverse local conditions on the beam propagation in the linear regime of low probing beam 
power of several pW, and compare our experimental results with numerical simulations. In 
order to numerically model light propagation along the z direction of any photonic lattice, 
we consider the paraxial wave equation {2i kd z + A^ — kn^r^ (^$ sc )} A(f) — 0 for a slowly 
varying electric field amplitude A(r) with wave number k. In this equation, d Xi is the spatial 
derivation in X{ direction, A^ the transverse Laplacian perpendicular to z, n e the unmodulated 
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Figure 3. Light propagation in aperiodic Fibonacci photonic lattices. Intensity distributions at exit 
face of crystal experimentally observed (first row) and numerically calculated (second row) for input 
probe beam size wq = 14 pm. Columns correspond to input beam positions 1 to 5, as shown in (k). 


refractive index for extraordinary polarization, 7*33 the corresponding electro-optical coefficient, 
and <F SC denotes a light potential caused by the optically induced internal electric field that 
in our case holds the Fibonacci structure. In order to simulate light propagation, we choose 
parameters that match experimental conditions and use a split-step method to evaluate the 
wave equation. 

For beam propagation studies, a Gaussian probe beam of wq — 14 pm beam waist is launched 
at different input positions of a Fibonacci photonic lattice. Figure [3] shows five selected output 
results of light propagation through a lattice of 20 mm propagation length. These different 
input positions on and between lattice sites are marked from number 1 to 5 in Fig. 0 k) and 
relate to output distributions shown in the first to fifth column in Fig. [3] as indicated. Top-row 
images present typical output distributions experimentally observed at the exit face of the 
crystal [Fig. [3](a— e)] while the bottom row represents the corresponding distributions obtained 
numerically [Fig. [3](f— h)]. 

We notice a very good agreement between experimentally obtained and numerically simu¬ 
lated results. Moreover, a pronounced heterogeneity of output profiles for different launching 
positions can be found indicating that the appearance of excited modes is highly diverse. 
Naturally, the separation between incident and neighboring lattice sites has a very strong in¬ 
fluence on the propagation process as the coupling coefficients vary inversely to the distance [7 . 
Thereby, light is subject to be guided along waveguide lattice sites accompanied by coupling 
between adjacent waveguides, as well as to ballistic propagation for predominantly high spatial 
frequencies encountering low spectral amplitudes [cf. Fig. |T](e)] . 


4 Comparison with beam propagation in periodic lattices 

To substantiate the localization character of Fibonacci structures, we compare this kind of 
deterministic aperiodic lattice with the case of periodic lattice, both of 20 mm length. Last 
named lattices are arranged periodically with as many lattice sites as for the aperiodic case, 
using again the sequential Bessel beam writing technique to preserve comparability. For this 
geometry, the lattice period corresponds to the effective distance of the Fibonacci lattice, 
d — 32 pm. The experimental effective intensity, the plane-wave probing output, and the 
focused Gaussian beam output are shown in Fig. |4j[a— c), respectively. Again, probing the 
structure with a plane wave does not resolve the light potential in very detail as can be found 
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in Fig. Qb), though, the image indicates qualitatively the area of modulated refractive index. 
A Gaussian beam of wo = 14 pm coupled into the central lattice point reveals the typical 
discrete diffraction at the output face as expected in regular lattices (7j. 

We further numerically study the beam propagation in Fibonacci and regular photonic lat¬ 
tices along the direction of propagation (z E [Omm, 40 mm]) with respect to the development 
of the probe beam width, having the experimental case of z = 20 mm included. In such kind of 
deterministic aperiodic as well as in disordered systems, it is common to perform analysis by 
averaging over different input beam positions [cf. Fig. p[k)] in order to diminish effects of the 
local environment given by the surrounding conditions |28,30,32|. For quantitative analysis, 
we use a relevant measure to characterize the degree of light confinement, the effective beam 
width cj e ff(z ) = P[z ) -1//2 , where 


= Jl-E(x,y,z)| 4 tfcdy 
(f \E{x,y,z)\ 2 dxdyf 

is the inverse participation ratio |28|. To extract an averaged effective beam width for the 
Fibonacci structure, we took the arithmetic mean of beam widths for 36 different on-site input 
positions. 

Figure [5] compares our numerical results of beam propagation in periodic lattices and ape¬ 
riodic Fibonacci photonic lattices. In Fig. [5](a) , averaged effective beam widths along the 
propagation distance z for Fibonacci and periodic lattice are presented for wq = 14 pm, which 
corresponds to ou e fi(z = 0) « 25 pm. The underlying numerical calculations consider a 15 x 15 
lattice site arrangement in order to avoid surface influences. The results in Fig. |5^a) indicate 
that c broadens slower during propagation in Fibonacci lattices than under periodic condi¬ 
tions. While all effective beam widths are almost equal for the first 10 mm, they separate for 
longer propagation distances since light starts to couple to next-nearest neighbors with respect 
to the input waveguide. At this point, aperiodicity comes into play at the earliest. 

Figure [5^b) presents the output after 40 mm of propagation in a periodic lattice. A typical 
signature of discrete diffraction is prominent here. The average of output intensity distributions 
for 36 different incident positions in Fibonacci lattice is presented in Fig. (5|c). Compared 
with the periodic case, localization in Fibonacci lattices is evident, as highest intensity is 
predominantly located around the input center. Such localization is weaker than expected 
in random media |l6 and diffraction is essentially suppressed in comparison with periodic 
lattices. 

Again, the remarkably heterogeneous behavior for different input positions in the Fibonacci 
lattice can be found in Fig. [5](a) where particular arrangements promote beam widths that are 



Figure 4. Experimental images for periodic photonic square lattice, (a) Effective intensity, (b) plane- 
wave, and (c) single-site probing output. 
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Figure 5. Comparison between strength of localization in periodic and Fibonacci lattice, (a) Beam 
widths cj e ff versus propagation distance z, (orange) mean beam width and (faint gray) single events in 
Fibonacci lattice; (dashed black) beam width in periodic lattice. Dashed vertical line at 2 : = 20 mm 
indicates experimental propagation distance. Intensity distributions at 2 = 40 mm in (b) periodic lattice, 
(c) averaged intensity, (d) maximum, and (e) minimum effective beam width in Fibonacci lattice, each 
for Wq = 14 pm. 


even larger than for the periodic case, especially around z = 30 mm. In contrast, configurations 
holding effective beam widths considerably below the mean value naturally result from most 
localized modes. To underline the diversity of localized modes occurring, we depict in Figures 
Id. e) conditions of maximal and minimal together confirming the relevance of the local 
waveguide arrangement conditions on light transport. 

Obviously, transverse light transport is slowed down in average due to aperiodic conditions 
accounting for localized modes whose appearance indirectly indicate the presence of band 
gaps [221. Periodic structures, in contrast, bring forward extended propagation modes as more 
light is carried away from the input position due to discrete diffraction. 


5 Conclusions 

To conclude, we have observed light localization in optically induced two-dimensional Fibonacci 
photonic lattices that arises from the aperiodic property of the waveguide arrangement. By 
experimentally and numerically analyzing linear propagation characteristics for various incident 
positions, we observed enhanced localization in Fibonacci photonic lattices compared with 
discrete-diffraction driven light transport in equivalent periodic realizations. We identify this 
localization with the occurrence of highly localized modes due to aperiodic structural conditions 
with diminished order, indirectly indicating the existence of band gaps. 




















A very good agreement between experimental and numerical results allowed to additionally 
analyze the development during propagation of a wave package sent to corresponding lattices. 
These numerical results consolidated our experimental observations and, moreover, gave deeper 
insights into propagation dynamics. We further found that selecting the input position in such a 
deterministic aperiodic structure is highly crucial to the beam broadening during propagation. 
Both, broader and smaller effective beam widths could be found among different scenarios 
compared to propagation in periodic waveguide arrangements. Evidently, the shape of excited 
localized modes and thus the strength of localization are local rather than global properties in 
these kind of structures that have been subject of our investigations. 

We are convinced that our results can be generalized to other kinds of aperiodic refractive 
index lattices, using the presented ideas and methods. 
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